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ABSTRACT 

We searched for infrared counterparts to the cluster of X-ray point sources discovered by Chandra 
in the Galactic Center Region (GCR). While the sources could be white dwarfs, neutron stars, or 
black holes accreting from stellar companions, their X-ray properties are consistent with magnetic 
Cataclysmic Variables, or High Mass X-ray Binaries (HMXB) at low accretion-rates. A direct way 
to decide between these possibilities and hence between alternative formation scenarios is to measure 
or constrain the luminosity distribution of the companions. Using infrared (J, H, K, Br7) imaging, 
we searched for counterparts corresponding to typical HMXB secondaries: spectral type BOV with 
K<15 at the GCR. We found no significant excess of bright stars in Chandra error circles, indicating 
that HMXBs are not the dominant X-ray source population, and account for fewer than 10% of the 
hardest X-ray sources. 

Subject headings: Galaxy: center, X-rays: binaries 



1. INTRODUCTION 

With the high spatial resolution of Chandra and 
ground-based infrared (IR) instrumentation, the Galac- 
tic Center, heavily obscured by dust, is now accessible to 
a wide range of astrophysical investigations. A series of 
Chandra observations of the GCR (~10' region around 
SgrA*), totaling; 67 ksec, have recently been analyzed by 
IMuno et all l |2003D Thereafter. M03), revealing 2357 point 
sources. The 2-8 keV X-ray luminosity (10 31-33 erg s _1 ) 
and spectral index (F < 1) of the majority of the sources 
are inconsistent with normal stars, active binaries, or 
young stellar objects, none of which are normally per- 
sistently luminous enough at >2.5 keV to produce the 
observed fluxes in the presence of the high absorption 
column (Njf ~ 10 23 cm -2 ). The only known sources ex- 
pected to be detected given the absorption are compact 
objects accreting from binary companions. The luminos- 
ity range of the GCR sources covers magnetic CVs and 
non-Roche-lobe filling HMXBs accreti ng in a low mass- 
trans fer regime, e.g. X Per, A0535+26 ( Ne gueruela et all 
2000) . Quiescent low mass X-ray binaries (LMXBs) with 
neutron star (NS) or black hole (BH) primaries are possi- 
ble contributors but are generally too faint (BHs) or too 
soft (NSs). HMXBs are young objects (few My), while 
CVs are slow to form and long-lived. The two possibili- 
ties point to different evolutionary phases of the GCR. 

If HMXBs comprise a significant fraction of the GCR 
sources, they must be the products of recent, massive- 
star formation. Frequent star-burst events in the GCR 
are expected on ^20 My scales based upon the insta- 
bility of large gas clouds orbiting within 200pc of the 
nucleus ijStark et al.l 1200^1 and the presence of several 
young clusters of massive stars. The HMX B hypothe- 
sis wa s proposed and discussed in detail bv iPfahl et al.l 
(200l. 

iGrindlavl l)1985(l suggested that the Galactic Bulge con- 
tains the debris of tidally-disrupted globular clusters. 
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Were the GCR sources dominated by magnetic CVs, 
globular clusters would then contribute because CVs are 
over-produced per unit mass in globulars relative to the 
G alactic plane. 

Morris ( 1993) proposed a third scenario in which com- 
pact stellar remnants form in the Galactic bulge and sink 
toward the center of the dense mass distribution via grav- 
itational interactions. The presence of a super-massive 
black hole (BH) at the center accelerates the process, 
lea ding to a cent r al den sity cusp. This scenario fits with 
the IMuno et all (|2005l) (M05) discovery of a centrally 
concentrated overabundance of X-ray transients, as well 
as the overall 1/8 source distribution found by M03. 

All three channels may be active now or have been in 
the past, so by measuring relative numbers of the various 
end-products, we hope to piece together the composition 
and history of the GCR. In this Letter we present prelim- 
inary results of our deep IR survey targeted at measuring 
the HMXB fraction. The occurrence of IR counterparts 
consistent with highly reddened massive stars and gi- 
ants will constrain the relative abundance of HMXBs and 
CVs. This work is part of ou r Chandra Multi- wav elength 
Plane Survey (ChaMPlane: iGrindlav et"aLll2005H which 
is an effort to map the space density of all types of low- 
luminosity accretion sources in the Galaxy. 

2. OBSERVATIONS 

The GCR is heavily obscured with average Ay ^25 so 
counterparts can only be detected with IR imaging. The 
stellar density is also very high so that reaching a suffi- 
cient limiting magnitude to make meaningful constraints 
on population models demands excellent image quality. 
In June 2004 we observed the inner 10'xlO' of the GCR 
with the PANIC imager on the 6.5m Magellan (Baade) 
telescope at Las Campanas, Chile. PANIC has a 2'x2' 
field of view with 0.125" pixels. Under good seeing con- 
ditions (median FWHM ~ 0.5"), we obtained a mosaic of 
25 pointings in J, H, Ks, and Br7 filters using a 5-point 
dither pattern with 3 exposures at each position. Sky 
emission was removed using off-source exposures, and the 
images were reduced and stacked in IRAF. Photometry 
was performed on the stacked images using SExtractor 
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ijBertin fc Arnoutsl 11996(1 . Astrometric and photomet- 
ric calibration were performed using the 2MASS point- 
source catalog. 

We used a version of X-ray point source catalog of M03, 
kindly provided by M. Muno, listing J2000 co-ordinates 
with net count rate, flux, and source significance deter- 
mined in 3 energy bands SI (0.3-1. 2keV), S2 (1.2-2.5 
keV), and H (2.5-8 keV). The intermediate S2 band pro- 
vides some additional information, but for this Letter 
we restrict ourselves to discussion of GCR (H only) and 
foreground (SI) sources. 

3. EXPERIMENTAL DESIGN 

We set out to identify all potential IR counterparts 
to the M03 X-ray sources and apply statistical tests to 
measure their numbers relative to random positional co- 
incidence. Our scientific goals are to set stringent upper- 
limits on the number of counterparts and examine their 
photometric properties. 

The availability of X-ray spectral information and 
multi-band IR photometry permits a clean separation 
of sources situated in the GCR from those in the fore- 
ground. The GCR sources are characterized by X-ray 
spectra completely cut off at low energies due to ab- 
sorbing gas and dust. Of the sources in the M03 cat- 
alog, 80% are detected at >3<r in the H band only, while 
any source detected in the SI band must lie in the fore- 
ground. The spatial distributions of X-ray point sources 
and K-band stars follow similar radial profiles about the 
GCR. M03 demonstrated a surface density with 1/8 de- 
pendence with angular offset 9 from SgrA*. Similar re- 
sults have_j2eCT_rerjort£d^^ observations 
(e.g. ICatchnol e. Whit elock 7^^179901) . 

There are ~100 known Galactic HMXBs, of which 70% 
have Be star companions whose spect ral types lie in th e 
narrow range 09-B2 and peak at BOV l|Negueruelall l998'). 
We adopt the luminosity and color of a BOV to rep resent 
"typi cal" HMXB values (My=-4, V-K=-0.83; lAlTenl 
2000). Note that many non-Be HMXBs are more lumi- 
nous because they have super-giant companions. For the 
CV case, we adopt a K0V dwarf (M V =5 .9,V-K=1 96) a t 
the brighter end of the normal range l|Pattersodll998jh 
Adopting Ay =25 and D=8 kpc for the GCR, our tar- 
get HMXB should have an apparent magnitude K=14.4 
while the assumed CV would have K=21. We note the 
intermediate polar GK Per with bright sub-giant coun- 
terpart would have K^19. 

3.1. Photometric Completeness Curve 

Due to severe crowding in the GCR, we infer the star 
density from our observations by correcting for incom- 
pleteness i.e. the loss of faint stars due to their brighter 
neighbors. Completeness curves for our PANIC pho- 
tometry as functions of magnitude and distance from 
SgrA* were determined by simulation. For each of the 
25 K-band stacked images, a random sample of stars 
was added using IRAF mkstar with a Gaussian PSF 
matched to the seeing. The photometry was repeated 
and the artificial stars identified and compared to their 
input magnitudes. A master catalog of input and recov- 
ered stars was then assembled, from which completeness 
curves could be extracted for any region of the 10'xlO' 
mosaic. The 50% completeness limit (fraction of stars 
recovered within ±0.5 mag of their input magnitude) 
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Fig. 1. — Star density curves for PANIC K-band photometry, 
in 5 magnitude ranges as indicated by legend. Observed number- 
density is shown by smooth lines, completeness corrected shown by 
histogram. A radial profile of the form A+Bx0~ 1 is also shown. 
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Fig. 2. — Peak-Up test for hard and soft X-ray sources, using 
PANIC K-band photometry. Solid and dotted lines show projec- 
tions along the RA and Dec axes of the number of X-ray sources 
with one or more stars inside the 95% error circle. A prominent 
peak is seen at zero offset only for the soft sources. The random 
match level and its fluctuations are apparent beyond 1" offset. 



within the inner 1' radius occurs at K=14.5 and over 
the rest of the mosaic at K=15.4. Curves for star den- 
sity are presented in Figure ^ and completeness values 
for the full field reported in the final column of Table ^ 

3.2. Peak-up Test 

The IR and X-ray catalogs are repeatedly matched over 
a large grid of positional offsets using integer multiples 
of the PANIC pixel-size of 0.125". At each offset posi- 
tion, the number of X-ray sources with one or more stars 
lying in their error-circle is recorded in a series of mag- 
nitude bins. We used 95% c onfidence er r or-circ les deter- 
mined from the formula of iHong et alJ l|2005|) which is 
dependent on net counts and distance from the Chandra 
aimpoint. For the GCR sources, error radii are 0.3"- 
1". The presence of real counterparts is signaled by a 
peak at the offset co-ordinates corresponding to perfect 
alignment of the two catalogs, providing a boresight. To 
confirm that X-ray and IR catalogs are on the same as- 
trometric frame we centroided the soft-source peak, find- 
ing a residual offset of ARA=0.011(1)", ADec=0.092(l)" 
with FWHM=0.55". 

The number of random matches (N ran ) is calculated 
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as the mean of all trials at shifts greater than 2", i.e. far 
away from correct alignment, and the standard deviation 
(a) of these values is used to characterize the uncertainty. 
We then take the (peak) number of matches (N m ), and 
subtract the random level to obtain the excess-above- 
random (N exc ). We claim detection of a population of 
counterparts above that predicted to occur by chance if 
N exc exceeds 3 x a. The upper-limit on the total num- 
ber of true counterparts was determined by monte-carlo 
simulation. At each iteration the X-ray catalog was off- 
set by randomly generated shifts in X and Y such that 
the absolute shift was between 2" and 5". A randomly 
selected sample of N trU e X-ray sources were flagged as 
"true counterparts" and the catalogs matched as previ- 
ously described. Since all of the X/IR matches were by 
definition random (due to the offsets) we were able to 
count the total number of matches N S i m in each mag- 
nitude bin, correcting for duplicates (sources flagged as 
simulated matches which also happened to have a ran- 
dom IR counterpart) such that N s i m — N ran + N trU e 
- N r an. and. true- Guided by the values of a from the 
peakup test, we ran simulations with A true =5-40 with 
100 iterations at each value. For each value of N trU e we 
recorded the largest value of N sim reached in 90% of tri- 
als, and fitted the results with a straight line relationship. 
Thus the upper-limit on the number of true matches cor- 
responds to the value of N trU e when N s i m = N m , these 
values are given in Table as (7go . 

4. CONSTRAINTS ON IR/X-RAY COUNTERPARTS 

We performed the peak-up test first on the entire field, 
and then independently on a series of concentric annu- 
lar regions centered on SgrA*. This allows us to limit 
the effects of crowding in the IR, and to follow the cusp 
distribution of hard X-ray sources. 

For the full field we found no significant excess in the 
number of IR matches for hard X-ray sources, irrespec- 
tive of magnitude, while for soft sources we found a 
highly significant (~30cr) counterpart population dom- 
inated by bright (K~ll) stars. This result is illustrated 
visually by Figure Eland in more detail in Tabled There 
is a marginal signal (2.9a) for hard sources versus the en- 
tire PANIC catalog, which could not be associated with 
any specific magnitude range. 

We therefore conclude the GCR counterpart popula- 
tion is dominated by stars fainter than our confusion 
limit, although some fraction of the hard matches are 
undoubtedly real. The 90% confidence upper limit for 
the full catalog indicates there are up to 87 hard counter- 
parts although we cannot at this stage identify them from 
among the 399 total matches. By placing completeness- 
corrected upper-limits on the number of counterparts as 
a function of K magnitude, we can constrain the propor- 
tion of HMXBs, and hence place a lower-limit on the pop- 
ulation of CVs and LMXBs. The final column of Tabled 
gives the cumulative percentage obtained by adding the 
upper-limits for each magnitude range and dividing by 
the number of X-ray sources. We find that a maximum of 
5.4% (at 90% confidence) of the sources could have coun- 
terparts brighter than K=15, given average GCR extinc- 
tion. The presence of a significant HMXB population 
would appear between K= 14-16, and the fact that we 
do detect soft-source counterparts for K=14-15 (at 6.8a) 
demonstrates that we are not merely being defeated by 
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Peakup Test on the full 10'xlO' field. 
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In the ~ 10 xlO PANIC survey area we consider: 
stars, 1453 hard X-ray sources, 105 soft X-ray 



sources. 2MASS point-source catalog used for K<12. N 7n = 
Number of X-ray sources with one or more matches. N ran = random 
matches, a uncertainty in number of random matches. N exc = 
number of "excess" matches above random level. Sig significance 
of excess. C Completeness fraction for K-band stars. Ugo= 
Completeness corrected upper-limit on the number of counterparts 
at 90% significance. U c = Cumulative Uqq brightward of K2 
expressed as 7, of number of X-ray sources. *No completeness 
correction applied for unbounded mag ranges. 



crowding in this range. 

For annular bins 0-1', T-2', 2'-3', 3'-4', 4'-5', our re- 
sults constrain the HMXB population as a function of 
distance from SgrA*. The effect of the nuclear cluster is 
apparent in the inner 1', where a large number of massive 
stars raise the upper-limits such that K<15 HMXB can 
potentially account for 9% of the 110 hard X-ray sources 
present (at 90% confidence). The limit is 5-8% for the 
remaining annuli. Overall the annular results are less 
constraining than the full field results due to the smaller 
sample sizes involved. A single annular bin (2'-3') pro- 
duced a peak at 3.8<r significance, possibly corresponding 
to a drop in the confusion limit sufficient to allow coun- 
terparts to be isolated. Such an effect could be due to 
seeing, anomalous extinction or an unrelated star cluster. 

We can account for 84% of the 105 soft sources down 
to K<15, and essentially all to K<16 where crowding 
dominates. The majority of soft sources will be dMe and 
other foreground coronal stars. For example at spectral 
type MOV our 50% completeness limit of K^15.4 implies 
a distance of just lkpc, indeed 86% of the SI sources 
have count e rpart s from our optical ChaMPlane survey 
iZhao et alJ l|200l . 

5. PHOTOMETRY OF POTENTIAL COUNTERPARTS 

Color-magnitude diagrams (CMD) were constructed to 
compare the H-K vs K distribution for stars lying in X- 
ray error-circles with that of the field population. Extinc- 
tion due to dust lying between us and the Galactic Center 
makes the GCR stars appear systematically redder than 
foreground stars, largely overcoming any degeneracy be- 
tween intrinsic brightness and distance. In figure|2|we see 
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Fig. 3. — Color-Magnitude diagram comparing the IR coun- 
terparts of hard and soft X-ray sources with the field population. 
Total numbers in each class given in histogram at right. (Satura- 
tion begins at K<11.5) 



the soft counterparts are mainly un-reddened stars, while 
the hard counterparts cluster around H-K~1.6 with a 
large scatter. Intrinsic H- K colors do not exceed ~+0.3 
for any stars {Allen 2000). Soft-source counterparts are 
also brighter than average field stars, the K-S probabil- 
ity that they are the same is 2.6xl0~ 9 . This result is 
consistent with the X-ray selection and peak-up test, the 
soft sources are shown to lie in the foreground, while the 
hard sources are dominated by random matches to stars 
in the GCR. 

Matches were found for 3 out of t he 8 periodic X-ray 
sources found bv lMuno etafl l)2004|) (7 in PANIC field, 
5 observed in Ks, see Fig. 0). One of these (CXOGCR 
J174534. 5-290201) is heavily reddened, while the others 
have slightly lower H-K values than the main concen- 
tration of stars in the CMD. The colors of the candi- 
date stars do not agree with the pred i cted extinction 
from X-ray spectral fits of lMuno etafl 1)2004) . The X- 
ray derived N# of >10 23 cm -2 requires the stars to be 
extremely bright to explain the observed K magnitudes 
and predicts H-K colors 2 mag redder than what is ob- 
served. These observations may be reconciled if most of 
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the X-ray-derived N# arises in the binary, in an accretion 
column or disk obscuring the X-ray source. A column 
of a few 10 22 cm -2 is typical for HMXB pulsars, but is 
insufficient to make up the difference. It is also possible 
that these matches are random and the real counterparts 
are much fainter. No match was found for the 7 bright 
transients detected by M05. 

A preliminary search for Br7 emission objects among 
X-ray IR counterparts yielded 11 candidates with 
S/N>5, of which 3 are foreground objects. Br7~K is 
severely restricted by the small amplitude of the expected 
signal (—0.1 mag for EW(Br7)~10-20A) and the effects 
of crowding. 

6. CONCLUSIONS 

The highly absorbed GCR X-ray sources appear to be 
dominated by a population of accreting binaries in which 
the mass-donor is less massive than in a typical HMXB. 
This conclusion follows from the lack of a significant ex- 
cess of stars brighter than K<15 in the error circles of 
hard X-ray sources, which is the apparent magnitude of 
a B0V star at the distance and mean Av of the GCR. By 
applying our peak-up test as a function of magnitude and 
distance from SgrA*, and correcting for completeness, we 
set 90% confidence upper limits on the number of hard 
X-ray sources that can have counterparts in a series of 
ranges between K=9-16. After accumulating these limits 
for all counterparts brighter than K=15, such stars can 
account for at most 79 out of 1453 hard X-ray sources, 
or 5.4%. Pushing one magnitude fainter to include Be 
HMXBs as late as B2, we find an upper-limit of 103 or 
7.1%. Within 1' of SgrA* these upper-limits become 9% 
and 15% due to the high density of very luminous stars 
in the nuclear cluster. We note that large variations in 
A„ and/or star clustering could modify these numbers. 

Our constraints imply non-HMXBs account for more 
than 90% of the M03 X-ray sources (at least beyond 1') 
favoring the leading alternative hypothesis of magnetic 
CVs. High angular resolution spectroscopy will be essen- 
tial to identify the true counterparts which likely exhibit 
B7 or He emission lines indicative of an accretion disk. 
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